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ABSTRACT 



The discovery of the first naturally occurring high gain hydrogen recombination line (HRL) maser, 
in the millimeter and submillimeter spectrum of the emission line star MWC349, requires an expansion 
of current paradigms about HRLs. In this paper we re-examine in general the physics of non-LTE 
populations in rocombining hydrogen and specify the conditions necessary for high-gain masing and 
lasing in HRLs. To do so we use the extensive new results on hydrogen level populations produced by 
Storey and Hummer (1995), and our calculations for the net (that is, line plus continuum) absorption 
coefficient for the hydrogen, and we present results for the a- and /3-lines whose principal quantum 
numbers n are between 5 and 100, for gas whose electron number density 3 < logiVg (cm~^) < 11, at two 
electron temperatures, Tg = 5, 000 and 10,000 K. We show that the unsaturated maser gain in an HRL 
is a sharp function of iVg, and thus to achieve high-gain masing each line requires a sufficiently extended 
region over which the density is rather closely specified. 

Saturation of masing recombination lines is a critical consideration. We derive a simple equation for 
estimating the degree of saturation from the observed flux density and the intcrferomctric and/or model 
information about the amplification path length, avoiding the vague issue of the solid angle of masing. We 
also present a qualitative way to approach the effects of saturation on adjacent emission lines, although 
the detailed modeling is highly case-specific. 

We draw attention to another non-LTE phenomenon active in hydrogen: the overcooling of popula- 
tions. This occurs for HRLs with n > 20, in gas where Ng < 10^ cm~^. Observationally the HRL overcooling 
might manifest itself as an anomalously weak emission recombination line, or as a "dasar," that is, an 
anomalously strong absorption line. In the simplest case of a homogeneous HII region the absorption 
can be observed on the proper free- free continuum of the region, if some conditions for the line or/and 
continuum optical depths are satisfied. 

We briefly discuss the prospects of detecting hydrogen masers, lasers and dasars in several classes of 
galactic and extragalactic objects, including compact HII regions. Be or Wolf-Rayet stars, and AGNs. 

Subject headings: masers — hydrogen: recombination lines — stars: WR stars — stars: Be stars — 
galaxies: AGN — stars: particular :MWC 349. 
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1. INTRODUCTION 



Masing spectral lines are an efficient observational probe in astrophysics: because of their strength 
and narrowness they often provide the best opportunity for detailed studies of kinematics and structure of 
the emitting source. Until recently only molecular astrophysical masers have been known. Yet calculations 
of the atomic hydrogen level populations in HII regions which have been done since the 1930's in order to 
interpret optical hydrogen recombination lines (HRLs) , showed that global population inversions probably 
exist across the Rydberg levels, thus raising in principle the possibility of lasing and masing in the HRLs. 
Numerical results directly implying inversion are found in papers like Gillie (1936), Baker and Menzel 
(1938), and, later, in Scaton's series of papers starting in 1959. 

In spite of that, and also in spite of the fact that Menzel (1937) pointed out (even before the discovery 
of laboratory masers and lasers) that under non-LTE conditions the integral of energy absorbed by a 
line may in principle become negative, the possibility of lasing in the HRLs attracted no attention for 
many years. It was perhaps due partly to the deeply rooted presumptions that the correction factor 
for stimulated emission, 1 — exp{—hi'/kT) , is negligible for optical lines because of their high frequency, 
and that HII regions arc always optically thin except perhaps in the Lyman HRLs, so that no serious 
amplification effects are possible anyway. Only relatively recently did these assumptions change. After 
modeling of optical HRLs in extragalactic objects suggested HII regions of unusually high densities (up to 
Nf, ~ lO^^cm""^) and/or large extensions, Krolik and McKcc (1978) considered the theoretical possibility 
of lasing in HRLs. Smith etal. (1979), analyzing their observations of an unusually intense infrared 
HRL in a dense HII region around the young stellar object BN in Orion, proposed that if the inversions 
persisted to the high densities expected in the source they could lead to laser amplification of this line. 

In the mid-1960's radio recombination lines were discovered. For them the correction factor for 
stimulated emission is much more important than for the optical and IR lines. Goldberg (1966) showed 
that even in HII regions having optical thickness |t| < 1 in both a radio recombination line and the 
adjacent free-free continuum, the inversion of populations can increase the line-to-continuum intensity 
ratio by as much as a factor of a few. 

The first high gain HRL maser was discovered in the millimeter recombination lines in the emission 
line star MWG349 (Martin-Pintado etal. 1989a). At first the maser emission was ascribed to the ionized 
outflow from the star (Martin-Pintado etal. 1989a,b), but later it was partly (Ponomarev etal. 1989), 
and then fully, attributed to the circum,stellar disk (Gordon 1992; Plancsas etal. 1992; Thum etal. 1992). 
Subsequent observational and interpretational papers enlarged on the list of observed masing lines in this 
source, and elaborated on the model in terms of a partly ionized masing Keplerian disk [Martin-Pintado 
etal. 1989b; Thum etal. 1994a,b; 1995; Gordon 1994; Ponomarev etal. 1994 ("Paper 1"); Strelnitski etal. 
1995a; 1995b ("Paper 3")]. 

With the discovery of real, high gain, HRL masers in MWG349, many basic theoretical questions 
arise about the physics of the pumping mechanism that maintains the inversion, the dependence of maser 
gain on the physical parameters (density, temperature, gas velocity) of the recombining gas, and about 
other possible non-LTE phenomena in HRLs. Directly related to these is the more practical question of 
discovering these non-LTE phenomena in other lines and in other sources. Though briefly touched upon 
in some papers (e.g.: Walmsley 1990, Ponomarev etal. 1991), these questions deserve to be studied on a 
more systematic basis. 

In this paper we readdress the theory of non-LTE phenomena in HRLs. Using the nomenclature 
for the non-LTE states of a quantum transition originally proposed by one of us (Strelnitski, 1983), 
along with the recent and extensive calculations of hydrogen level populations by Storey and Hummer 
(1995), we analyze this phenomena (Section 2). We demonstrate that the population calculations predict 
not only inversion and its weaker complement, "overheating" , both of which increase the intensities of 
emission lines, but also their opposite — "overcooling" of populations — which tends to produce an 
anomalously strong absorption in the line (the "dasar" effect). We briefly discuss from a physical and 
thermodynamical point of view how these non-LTE states are created in recombining hydrogen. 

We address in Section 3 the conditions necessary for producing high maser gain in an HRL, and the 
conditions (and observational criteria) of saturation . The "net" absorption coefficient which accounts for 
both the negative absorption (amplification) in the line and the positive free-free absorption at the line 
frequencies, proved to be a useful parameter (Ponomarev 1994). Our calculations of the net absorption 
coefficient allow us to derive the optimum electron density Ng for strong unsaturated masing in a- and 
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/3-lincs with 5 < n < 100 for gas with densities up to A'',. = 10^^ cm^'', and at two representative vahies 
of electron temperature, Tg = 5,000 and 10,000 K. The results are illustrated graphically and tabulated. 
Rough analytical approximations are provided for the lower n transitions. We derive a formula for 
estimating the degree of saturation in an observed masing HRL, if the line is formed at its optimum 
electron density. 

Section 4 is devoted to the phenomenon of overcooling and "dasing" in HRLs. We use the solution of 
the line-plus-continuum radiative transfer equation for a homogenous slab to show that direct evidence of 
overcooling of a line would be its appearance in absorption on the free-free continuum of the HIT region. 
We derive the minimal optical depths in the line and/or the continuum necessary for this to happen. 
This is not, however, the only way to detect overcooling: the increment of HRLs seen in emission can 
demonstrate it as well. 

Finally, in Section 5, we discuss the prospects of detecting HRL masers, lasers and/or dasers in 

galactic and extragalactic sources other than MWC349. MWC349 itself the only high-gain HRL 
maser detected so far — is addressed separately in a forthcoming paper (Strelnitski et al. 19956, "Paper 
3"). 



2. NON-EQUILIBRIUM STATES OF RYDBERG TRANSITIONS 
IN ATOMIC HYDROGEN 

2.1. Radiative Transfer 

Our analysis of non-equilibrium effects in HRLs is based on the standard solution of the radiative 
transfer equation for the simplest case of a steady-state, homogeneous, plane-parallel medium of finite 
optical thickness t^: 

I, = S,[l- exp(-T,)] . (2.1.1) 

Here is the outgoing radiation intensity at the frequency v and is the source function. No external 
background source is assumed. In the case of the HRLs, where both line and free-free continuum are 
formed together, it is convinient to represent the sorce function as = 'qvB„{Te) (Goldberg 1966), 
where 

- B,{T,) k', + fc;6i/3i2 ' ^ - ' ' 

Bv{Te) is the Planck function for the electron temperature T^; k'^ and k'^ are the continuum and the 
line absorption coefficients at LTE, both corrected for stimulated emission (the prime indicates LTE 
values); &„ is the MenzePs coefficient of departure of the level n population from its LTE vahie. /?i2 in 
Eq. (2.1.2) is formally the ratio of the actual (non-LTE) and the LTE correction factors for stimulated 
emission (Brocklehurst and Seaton 1972): 

1 - (62/61) eXp(-/lt/o/fcTe) (r, 1 nx 

1 — exp {—huo/kTe) 

where Uq is the rest frequency of the 1-2 transition and constants have their usual meanings. /?i2 is a 
measure of the deviation of the population ratio of two levels from its LTE value. 

The net optical depth at the line center for a homogeneous medium is (we omit the subindex vq): 

Tnet =Tc + Tl= k'^L + k'l bi(3i2L = (fc^ + k) L = knetL . (2.1.4) 

where fc; = k[ 6i/?i2 is the actual (non-equilibrium) line absorption coefficient and fc„et is the "net" 
(continuum plus line) absorption coefficient. 
Since 

= B^Te) [1 - exp(-re)] , (2.1.5) 
the relative line intensity at the line center is given by: 

r (^o) ^ '^""i 7 'f'^ ^ V \ '"""^r^i - 1 = / : '"'^ - 1 • (2.1.6) 

Ic{i^o) 1 - exp(-Tc) 1 - exp(-Tc) 
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All r's and r] here and following are at the central line frequency vq, with the subindex omitted. 



2.2. Non- Equilibrium States of a Quantum Transition 

There is another commonly used parameter, besides /3i2, to describe deviations of a quantum tran- 
sition from equilibrium: the excitation temperature T^, which is defined by the Boltzmann-like formula 
for the population ratio of the two levels: 

= = exp i-huo/kT,) , (2.2.1) 

or 

- - (2.2.2) 



ln(ni/n2) 

In this subsection and in Fig. 1 we denote by rii the population density per one degeneracy sublevel: 
m = Ni/gi, with Qi being the degeneracy of the level i; this should not lead to a confusion with the 
principal quantum number for which the script n (without subindex) is used throughout the paper. 
Subindex 1 denotes the lower and 2 the upper level of the transition. 
Since 

n2 _ 62 n'2 
ni bi n'l ' 

where primed symbols denote again the LTE values, we have the following relation between (3i2 and T^: 

_ 1 - exp{-huo/kT^) 



For low frequency transitions {hi/o/k <gC Tg ), and when also \Tx\ t 



/3i2 « ^ . (2.2.4) 

The behavior of and /9i2, as functions of nijni. is illustrated in Fig. 1 for the case of T,, = 
10, 000 K and hv^/k = 100 K. When the population ratio ni/n2 decreases from 00 (all the atoms in 
the lower level) to (all the atoms in the higher level) , first increases from to 00 , then undergoes 
a jump to —00 at ui = n2, and then increases to 0. The 0i2 parameter decreases monotonicly from 
pmax _ 1/(1 — exp(— /ir/o/AiTe)) to —00, passing through at ni = n2- 

The classification of the excitation states of a transition is most natural in terms of temperatures 
(Strelnitski 1983), see Fig. 1. = T,. is obviously the state of thermalization, and /3i2 = 1 in this 
state. When < < +00, we refer to the condition as one of overheating of the transition 1-2; it 
corresponds to < /?i2 < 1. The case <0 (with /3i2 also being negative here) is the case of population 
inversion. Finally, when > T!,. > (corresponding to /3i2 > 1), we call this condition overcooling 
of the transition 1-2. The choice of this nomenclature will become clearer as we analyze the observable 
effects of the corresponding non- LTE states below (see also Strelnitski 1983). 



2.3. Non-Equilibrium Populations in Hydrogen 

Storey and Hummer (1995) have recently completed a detailed set of calculations of the hydrogen 
level populations for levels n < 50, and presented the results as machine-readable files. For Menzel's 
Case B (the Lyman lines are infinitely optically thick, all other lines are thin), calculations were made for 
2 < logiVe < 14 and for 2.7 < logTe < 5 . Figs. 2-5, plotted with the use of the Storey and Hummer's 
data, show that recombining hydrogen may acquire all the types of non-LTE population distributions 
just described overcooling, overheating and inversion. 

Fig. 2 presents the 6„ coefficients versus principal quantum number n, and as function of electron 
density, A^e, for = 10^ K. The major qualitative features of the 6„ behavior are: 



5 



(1) /;„, tends to one when n tends to oc. This thermahzation of the highest n levels is d\ie to the fact 
that their populations are dominated by three-body coUisional recombinations and coUisional ionizations, 
both processes involving the same Maxwellian gas with kinetic temperature Tg. Thus, a reservoir of 
quasi-thcrmalizcd levels is created at the highest n's. 

(2) hn decreases toward lower n, due to the ever increasing rate of spontaneous decay which is roughly 
proportional to . There is a diffuse separation between those levels which are quasi-thermalized and 
those which arc underpopulated and show a steep positive gradient dfe„/dn. This boundary progressively 
shifts toward lower n's, as A^e (and the effects of collisions) increases. 

(3) At the lowest n's the gradient d6„/dn tends to some high positive limiting value. Here populations 
are totally controlled by radiative processes, and the intrinsic trends of radiative probabilities across the 
levels determine this "radiative limit" in the slope of (compare this part of the 6„(n) curves with 
our bn calculation for the case A^e = 0, shown in Fig. 2 with the broken line). The passage to the radiative 
regime, at low n, occurs because, at a given density, the rate of transitions induced by collisions with 
electrons decreases with decreasing n approximately as n**^^, whereas the rate of spontaneous radiative 
decay increases as n~^. Thus the relative role of radiative processes, at a given density, increases with 
decreasing n very steeply, steeper than . 

(4) At low enough densities (/^e ~ 10^ cm~^ ), d6„/dn changes sign before passing to the radiative 
limit. As seen from Eq. (2.1.3), this means that (3\2 exceeds +1, and thus the corresponding transitions 
are overcooled. To our knowledge, this important fact has never been pointed out before, although it 
has been implicitly present in the results of earlier (n, /) population calculations. It is clearly seen, for 
example, in Fig. 2 of Hummer and Storey (1992). The change of sign of d6„/dn is due to the onset 
of partial "unbhirring" of the levels' angular momentum /-structure popultations at these values of n, 
a distribution which is blurred by proton collisions at higher n. We interpret it as follows. Unblurring 
begins when the radiative decay rate of an (n, /) sublevel becomes comparable with the collision rate with 
protons. Since spontaneous decay rates increase with decreasing /, the low I sublevels of a given n level 
are unblurred the first. At any given density, therefore, there is a transitional group of levels in which 
more and more Z-sublevels (beginning with the smaller I values) are unblurred as one looks to smaller n. 
Blurring of the /-sublevels enables high / atoms on a level n to decay faster, via lower / sublevels, than 
the A/ = ±1 selection rule normally allows. The more a level is unblurred the higher is the / which looses 
this ability to decay more rapidly. As a result, the lower the principal quantum number n the more its 
bni factors reflect higher populations in its higher / siiblevels. This effect is seen in Fig. 3, plotted with 
the numerical data of Storey and Hummer (1995) for = 10, 000 K and N,, = lO'^cm"'^. The effective 
bn factor for a level n is derived from the statistically weighted, calculated 6„; factors: 



Since high / sublevels have higher statistical weights (2Z + 1), the increase of their bni factors is more 
significant for the weighted value of 6„ (2.3.1) than the dc;c;rease of the bj-a factors of low I sublevels. 
As a result, unblurring increases the effective bn value of the whole n-level. And because the degree of 
unblurring increases with decreasing n, a negative gradient d6„/dn may result for the partly unblurred 
group of levels, as is seen in Fig. 2. We note that those prc^vious calculations which neglected the angular 
momentum state population distributions (e.g. Walmsley 1990) show no evidence for "overcooling." 

These four features of the 6„ behavior determine the shape of the /3„,„+i vs. n curves for a-lines, 
shown in Figs. Aa.b for several values of A^e, and the shape of the Pn,n+i vs. A/g curves shown in Fig. 5 
for several a-transitions. Comparing Figs. 4a, 6 with Fig. 2 is instructive. 

(1) At the high values of n the positive slope of 6n vs. n results in an inversion of populations — 

I3n,n+1 < 0. 

(2) In a transitional set of n values 6„ levels off and then changes slope due to the unblurring of I 
states at progressively smaller n values; here the inversion first diminishes, and then disappears. If the 
density is low enough (Ne t cm^"^ ), the imblurring of I states first reduces the inversion to overheating 
(0 < (in,n+i < 1 )) then leads to thermahzation {Pn,n+i = 1 )> and finally — to overcooling (/3„,„+i > 1). 

(3) The final transition of 6„ to the "radiative limit," at still lower n, stops the trend towards 
overcooling and brings (in,n+i down to the slightly overheated values typical for the "radiative limit" at 
low n. 




(2.3.1) 
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At low densities, therefore, the /3n,n+i (^) vs. n func;tion has two extrerna a minimum due to 
inversion at high values of n , and a maximum due to overcooling at lower n's. At densities Nf. ^ 10^ cm~'^ 
there is only the inversion minimum because ion collisions at these high densities blur the /-structure 
down to very low n's, and thereby remove the effective cause of overcooling. 

Fig. 5 illustrates the behavior of the /3n,n+i factor as a function of A^e for some selected a-transitions. 
Qualitatively, the (3 (Ne) functions are similar to the /3 (n) functions: high n lines (higher than n w 20) 
have two extrema — a minimum due to inversion at high densities and a maximum due to overcooling at 
low densities. Lower n lines show only inversion, and this at increasingly higher densities (and diminishing 
amplitude) for lines of decreasing n. 



2.4- Thermodynamic Considerations 

The creation of stable conditions for inversion, overheating, or overcooling requires cycles of pop- 
ulation transfer which, from the thermodynamic point of view, operate as "heat pumps" (Strelnitski 
1983). The system of population transfer cycles in hydrogen is rather complicated but we can infer some 
general principles. The two ultimate high temperature and low temperature reservoirs that interact with 
the system of hydrogen levels and provide free energy to create and maintain non-equilibrium popula- 
tions arc, respectively, the "hot" ionizing UV radiation field and the "cold" radiation field at the lower 
frequency recombination lines. Note that this difference in the high-frequency and the low-frequency ra- 
diation temperatures arises not because the central star is hot: a diluted Planck spectrum with any color 
temperature possesses this gradient of radiative temperatures (Strelnitski 1983, Sobolev etal. 1985). In 
such a radiation field the net cycle of quantum transitions is always direct (distributive) — high energy 
photons are split up to lower energy photons — this constitutes the content of the Rosseland theorem 
(Rosseland 1926; sec also Strelnitski 1983, and Sobolev etal. 1985). This cycle of fluorescence makes 
up the grand radiative-radiative (RR) pumping cycle in the recombining hydrogen in HII regions, and 
tends to produce a global population inversion due to the steep increase of the sink (spontaneous decay) 
rate with the decreasing n. However, besides the ultimate source and sink radiative reservoirs, there are 
other, intermediate, reservoirs involved with the process: the recombination transition levels themselves, 
with their excitation temperatures, and the Maxwellized gases of electrons and protons with their kinetic 
tcmpcratiircs T,,. The heat exchange between these reservoirs provides secondary cycles of population 
transfer. Each of the two collisional reservoirs, free electrons and protons, have different degrees of effec- 
tiveness in their interaction with the different hydrogen levels and sublevels, and because of this difference 
they modify the pattern of the primary (RR) cycles and can even change, locally, the sign of the cycles, 
producing overcooling. 

We note that only one net link — the source link or the sink link — can be collisional in these 
secondary cycles; the other link must be radiative. In other words: pumping cycles in hydrogen can be 
of the radiative-radiative (RR), collisional-radiative (CR), or radiative-collisional (RC) type, but they 
can not be of the collisional- collisional (CC) type because only one temperature, Tg, is associated with 
collisions in HII regions. 

3. HYDROGEN MASERS 
3.1. Maser Gain 

The overheating and inversion of populations occur when (3i2 is smaller than its LTE value P'12 — 1- 
It is seen from Eqs. (2.1.2), (2.1.4) and (2.1.6) that 77 is then higher than its LTE value 77' = 1, T„et is 
lower than its LTE value, and in either case it increases the relative intensity of the line. Goldberg (1966) 
demonstrated that in this situation the relative intensity of the line r can increase up to a factor of a few 
even when \ti\ and Tc are < 1. 

However, to produce a maser amplification » 1, population inversion (3i2 < alone is not suflicient. 
The net optical depth Tnet must also be negative and its absolute value must be considerably higher than 
unity. When Tnet < 0) both 77 and the numerator in Eq. (2.1.6) are negative, and Eq. (2.1.6) becomes: 

^ ~ ~/ V *^^P \'^net\ ■ (3.1.1) 

1 - exp(-Tc) 
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Thus the relative intensity of a recombination Une may become very high, if the exponent |T77,ef; I (the 
"gain" of the maser) is » 1. 

The maser gain is determined by the length of the amplification path, L, and the value of knet = 
k'^ + k'lbiPu [sec Eq. (2.1.4)]. A negative value of knet can only derive from a negative /3i2, because all 
the other quantities determining it are essentially positive. Since k'^, in particular, is positive, free-free 
absorption always reduces the maser gain, with the amount of reduction depending on the frequency. 

If the density of an HII region is not very high (say, ^ 10^ cm~^), /3i2 for a given transition can attain 
high negative values at the bottom of its inversion minimum (cf. Fig. 5), thereby increasing significantly 
the absolute value of ki = fc;6i/?i2 as compared with its LTE value fc;. Yet, the density at which maser 
gain in a transition attains its maximum docs not coincide with the density of this maximum inversion. 
The reason is that both /3i2 and k'^ depend, and depend differently, on electron density N^, in particular 
with the latter having a quadratic dependence on Ng. Multiplication of /?i2 (A^e) by the steep (oc N^) 
function fc; (A^e) shifts the maximum of the product toward higher Ne- It also makes the range of densities 
with high \ki \ values much narrower than the range of densities with high \/3i2\ values. Furthermore, since 
ki and k'^ depend on A^e differently, accounting for free-free absorption shifts (downward) somewhat the 
value of the density where fc„et is minimum with respect to that for ki alone. 

In order to derive the optimum conditions for high-gain HRL masing we have performed extensive 
calculations of the knet i^e) dependence for a- and /3-lines with 5 < n < 100. Since the minimum value 
of knet for any line occurs at a density much higher than the density of unblurring the Z-structure, our 
program calculating 6„ coefficients neglects the ^-structure as unimportant. 

Several illustrative cases of the calculated knet i^e) dependence, for = 10,000 K, are shown in 
Figs. 6a f (a-lines) and 6g (/3-lines). In Fig. 7 the knet (Ne) dependence is given for a large interval of 
higher n. A logarithmic scale on the ordinate has been used to show all these lines on one graph. Only 
the locus of the maximal \knet\ values for individual lines and the envelope of the |fc„et| (A^e) plots are 
shown for Te = 5,000K. The knet (n) dependence for different densities, from A^e = 10^ to 10"^^ cm~^, 
and for the two temperatures is presented in Fig. 8. 

Comparison of the knet i^e) for the 36a line (Fig. 6b) with the 0n,n+i i^e) for this line (Fig.5; the 
intervals of high \ki\ and \knet\ are indicated by the bars) shows that whereas /?«,„+! is negative over a 
large density interval, with a maximum of inversion that peaks at A^e ~ 10^ cm~"^, the density interval of 
high \knet \ values for this line is much narrower, and peaks at much higher density — A^e ~ lO^'cm"^. 

Table 1 gives N™"'^, the electron number density at which |fc„et| attains its maximum, and the 
corresponding value of knet, for twenty a-lines and twenty /3-lines at Te =10,000 and 5,000 K. The 
dependence A^™"^ (n) is shown graphically in Fig. 9. For the lines with 40 (these are of particular 
interest because many of them have already been observed as strong masers in MWC349), we fit the 
following analytical approximation for TV™"^ (n): 

jynrax _ g q . -^q15 ^-5.66 (j.^ ^ ^q4 ^ (-3 -^ 2) 

^rnax ^ J J ,^Ql5^-5.75 (Tg = 5 • 10^ K) . (3.1.3) 

The A;^"^ (n) dependence is approximated in the same interval of n by: 

kZT « -2.3 • 10-3 n-« ° (Te = 10^ K) , (3.1.4) 
Cef « -3-4 • 10-2 n-^-^ (Te = 5 • 10^ K) . (3.1.5) 

The difference between these approximations and more exact model calcuations is under « 50% for 
5 < n < 40. The approximation (3.1.2) is drawn in Fig. 9 as a broken line. 

We point out that N™""^ for a given transition doesn't depend strongly on the assumed temperature. 
It is seen in Table 1 and in Eqs. (3.1.2) and (3.1.3) which are, in fact, almost identical. For rough 
estimates of iV™"^ in the whole temperature interval 5,000 10,000 K their "average," 

j^max _ g . ;^Qi5 ^-5.7 . ^q3, < j; < K) , (3.1.6) 

can therefore be used. 

In contrast, knei (N"'-""'') may vary significantly with temperature, especially at the higher n values. 
The temperature dependence of knet is more generally seen in Figs. 7 and 8. 

The narrowness of the knet i^e) peaks for individual lines demonstrated in Figure 6 implies that to 
reach high gain masing in any small group of recombination lines around some n an extended region at 
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the density close to the "optimum" one for this n is required. However, if a vahie of is fixed, rather 
than a value of n, the corresponding maximum value of the gain, |A:™j^| (A^e); is not that of the line 
attaining its maximum gain at this density: there are several lines, of smaller n, whose gain is higher; 
likewise, there are several lines, of higher n, whose gain is smaller. This can be seen in several graphs of 
Fig. 6, but especially clearly in Figs. 7 and 8: the locus of |fc„et| values for individual a- lines in Fig. 7 (the 
two broken lines, for two temperature values) passes below the locus of maximal \knet] values for given 
densities (the envelope of the individual plots). Note that both these loci are straight lines in the log |fc„et | 
- logiVe graph, displaying that corresponding dependencies are of the power-law character. Fig. 8 also 
shows that the lower the density, the flatter the maximum of the |fc„et| (n) function. Thus, in contrast to 
the case of a given line where a narrow density interval is required to attain the maximum possible gain, 
in the case of a given density (if the density is low enough, say, ^ 10'' cm~'^) many lines are allowed to 
mase with comparable gains. 

Li = \knet\~^ is the length corresponding to |r„et| = 1- Figures 6 8, Table 1 and Eqs. (3.1.4) (3.1.5) 
thus allow us to estimate how extended the medium of a given electron density should be along the 
line of sight to provide a maser gain > 1 in a corresponding HRL. We note, however, that all these 
give only an upper limit to a or a |fc„et| because the calculations were done assuming the local line 
width is only due to the thermal broadening, ignoring possible microturbulent broadening. If the latter 
is present and described by an equivalent temperature Tt^rb, the value of \ki\ decreases by a factor of 
[(Tg + Tturbj/Te]^^"^ ■ The relative decrease of \knet\ due to any microturbulent braodening (and the 
corresponding increase of Li) is even larger than that for \ki\, because of the contribution from k[. as seen 
in Eq. (2.1.4). The value of Li may additionally be increased by the presence of large scale motions in 
the source (expansion, contraction) which diminish the column densities of atoms at every radial velocity 
(frequency), and thereby increase the value of Li by a factor of approximately AV/Av, where is the 
dispersion of radial velocities due to the large scale motions and A?; is the local radial velocity dispersion 
due to the thermal motions and microturbulence. 

All the data for Fig. 6, Table 1, and Eqs. (3.1.2) - (3.1.5) were obtained assuming the media are 
optically thin in all recombination lines (except for the Lyman lines which, for the Menzel's Case B, are 
taken to be infinitely optically thick; we do not consider these lines here). Yet, since we deal here with 
high-gain masers, that is with gas which is supposed to be optically thick in masing lines, and since, in 
hydrogen, the masing and pumping lines are intermixed and have comparable values of absorption coeffi- 
cients, the assumption of optical thinness in all but masing lines is internally contradictory. Qualitatively, 
the effect of radiation trapping is similar in the "ordinary" (non-inverted) and in the inverted transition 
cases: the more the radiation is trapped (the larger the modulus of the line optical depth), the more the 
transitions induced by the trapped photons compete with spontaneous and collision-induced transitions 
in controlling the populations. The results of population calculations for a multilevel system with and 
without account for radiative trapping may, of course, differ drastically. It is impossible in practice to 
accoimt in a general way for radiation trapping, not so much because it is a mathematically complicated 
nonlinear problem, but rather because the pattern of optical depths is peculiar to each concrete source. 

Fortunately, in the important particular case of masing in a Keplerian disk (the case of MWC349 
and similar objects) the problem of the pumping radiation trapping can be finessed by supposing that 
the medium is optically thick in only one direction, that of the chord of maximum velocity coherence (see 
Fig. 1 in our Paper 1). In the directions perpendicular to the disk, and also in those directions within 
the disk which are far from that chord, the optical depth is smaller, and it can actually be smaller than 
unity, even if the chord optical depth is considerably larger than unity. Since these directions occupy, all 
together, much larger solid angle than the directions close to the chord, spontaneously radiated photons 
can escape from the medium almost freely, providing an efficient pumping for high-gain, optically thick 
masing into relatively small solid angle along the chord. There is one further fact providing confidence 
in this general approach: the successful modeling of the masing observed in the case of MWC349 (cf. 
Paper 3). 

If the modulus of the optical depth along the chord of maximum amplification is large enough, the 
rapid increase of maser intensity creates a special kind of radiation trapping problem — maser saturation. 

3.2. Saturation 
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Saturation ocx'urs when the growing mascr radiation begins to influenc;e the; populations of the 
transitions significantly. No general solution exists for the problem of maser saturation for the same reason 
why no general solution is possible for radiation trapping in a system with positive optical depth. Yet, 
as it will be seen below and in Paper 3, saturation may bo occuring in these systems, and understanding 
its role will play a crucial part in interpreting the observed properties of masing HRLs, and in predicting 
the observability of other lines. In this section we consider two basic questions concerning saturation in 
HRLs: (1) how saturation in a masing line affects the level populations (thus the gain) of the masing 
transition itself; and (2) how, in a multi-line HRL maser, saturation in one line affects the populations 
and the maser gains in adjacent lines. We also derive a criterion for suspecting when saturation has 
occurcd in a masing HRL. 

The first question has been intensively studied in the context of molecular masers. The condition 
of saturation is that the maser radiation density influences "noticeably" the population difference AA^2i 
between the upper (2) and lower (1) maser levels. Ideally one needs an analytical expression for AiV2i as 
a function of all the basic processes influencing the populations. In practice this goal is unattainable when 
the number of interacting levels is > 3 — 4. Several simplifying approaches have therefore been proposed 
to get around this problem (see the discussion in Strelnitski, 1993). In all these approaches the two 
relaxation processes, radiative (maser) and collisional transitions between the maser levels themselves, 
are accounted for almost exactly. It is much more difficult, however, to provide for a complex system a 
simple yet adequate approximation for all the pumping processes. 

In one approach, pumping is represented by four phenomenological coefficients for the upper and 
lower maser levels: the "income" coefficients A2 and Ai (cm~^s~^) and the "outlay" coefficients r2 and 
Fi (s""'^). Omitting subindex "12" where it cannot lead to misunderstanding, and taking for simplicity 
the degeneracies of the maser levels to be equal, the solution of the statistical equilibrium equations in 
this "A, F" -approach can be represented as follows (e.g.: Strelnitski, 1993): 



AN = ° = , (3.2.1) 

^ ^ r+C2i Js 



where 

A2 



ANo = (3.2.2) 

r + 621 

is the unsaturated population difference; F = (F2 ^ +T^^)~^ is the harmonic mean value of the decay 
rate of the two levels; B is the Einstein coefficient for induced radiative transitions between the maser 

levels {B = B12 = B21 since we assumed the degeneracies of the two levels to be equal); J is the actual 
radiation intensity averaged over the profile of the absorption coefficient and over directions; 

J, = {r + C2i)/B (3.2.3) 

is the so called "saturation intensity;" and C21 is the rate of collisional relaxation between the maser 
levels. To obtain this form of eqs. (3.2.1) and (3.2.2) we ignored the slight differences between C12 and 
C21. Strictly speaking, then, these equations are only valid for transitions with hu <C kTf. which, at 
Tg ~ 10^ K, comprise all hydrogen a-lines down to « 5a. 

Eq. (3.2.1) shows that maser radiation begins to influence the maser level populations when the rate 
of its interaction with the maser levels, BJ, becomes comparable to the higher of the two rates: the rate 
of collisional relaxation between the levels, C21, and the rate of the decay of the populations from the 
maser levels in the pumping processes, F. (It is worth noting that harmonic mean is always smaller than 
the smallest of the two averaged numbers, and equal to one half of each when they are equal.) The net 
effect of the mascr radiation is the transfer of population from the upper to the lower state. It is seen in 
Eq. (3.2.1) that in the limit J > the population difference AA^ 0, and thus {N2/N1) 1. Note, 
however, that in this limit, AA'' oc [see Eq. (3.2.1) with J ^ Jg]. This dependence of AA'' on J 
reduces the mascr amplification from an exponential [Eq. (3.1.1)] to a linear dependence. The effect is 
easily understood by noting that the radiative transfer equation for a sufficiently strong maser (where 
the spontaneous term is negligible when compared to the induced term) reduces to: 

^ = l^BANI. 
dz Ait 
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When AN is independent of the intensity / (unsaturated regime) . integration gives an exponential growth 
of / with z; when AA^ oc (saturation), I depends on z hnearly (see also general reviews of maser 
theory, e.g., Strelnitski, 1974; Reid and Moran, 1981; Elitzur, 1992). 

The major simplifying assumption of the "A, F" -approach is that X2 and Ai do not depend on the 
maser level populations. This leads eventually to the expression (3.2.1) which states that {AN/ANq) 
does not depend on A's but only on F's and C — a radical simplification. The fewer the number of the 
levels involved with the pumping, the worse the approximation, because in such a few-level system the 
non-maser level populations — the source of "income" for the maser levels — depend sensitively on the 
populations of the maser levels. 

In the case of HRLs pumped by recombinations, cascading and collisions, every level is fed by many 
other levels. If there is only one masing transition in the system, the change of its populations caused 
by saturation should not produce any serious distortion of the reservoir of pumping. Moreover, even if 
masing and saturation occur in several adjacent lines in the same region of space, this will hardly change 
the reservoir of pumping for individual transitions because saturation only re-distributes populations 
between two adjacent levels, without bringing populations in or out of the reservoir as a whole. Thus 
the "A, r" -approximation should give an adequate qualitative description of the AN{J) dependence for 
individual masing recombination lines, in spite of the fact that in hydrogen the reservoir of populations 
for pumping may heavily overlap with masing transitions. 

The hydrogen levels we consider here decay primarily through spontaneous and collision-induced 
transitions to other levels; thus, 

r^At + Ct, (3.2.4) 

where = (A^^ -|-^^^)~-^ and Cj = {C^^ +C^^)~^ are the harmonic mean values of the total Einstein 

A coefficient and of the total collision rate for the two maser levels respectively. Using Eq. (3.2.4) and 
the relation between the Einstein A and B coefficients, we can rewrite Eq. (3.2.3) as 

^ 2K' At + Ct + C2i 

where ^421 is the Einstein coefficient for the maser transition, z/q is the transition frequency and the 
constants have their usual meanings. 

To estimate the degree of saturation in an observed masing line, we have to compare an estimate 
of J in the source with Jg given by Eq. (3.2.5). A key and complex issue is that of the radiation's solid 
angle, but we are able to sidestep the matter using the following argument (cf., Strelnitski, 1984). On 
the one hand, unless the geometry of the source is very special, the solid angle of the maser beam should 
be close to 

J7 « ^ , (3.2.6) 

where a is the observed emitting surface area, and / is the amplification path length. On the other hand, 
the solid angle under which we see the source is 

(3.2.7) 

where D is the distance to the source. 

Thus, a cancels out in the equation connecting J with the observed flux density S: 

,.4 = 5"=5i^. (3.2.8) 

where / is the specific radiation intensity averaged over the profile of the absorption coeflacient. Prom 
Eqs. (3.2.5) and (3.2.8) we obtain the degree of maser saturation: 

Js Snhp^ At + Ct + C2iP' ^ ■ ■ ^ 

Eq. (3.2.9) gives an estimate of the degree of saturation in terms of the transition constants 
vo, A21, At, the model parameters Ct, C21 (determined by the model density of the emitting region), 
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the observed flux parameter S. and the "semi-observed" parameter / which, though not directly observ- 
able, can in principle be reliably estimated from direct interferometry of the source structure in the plane 
of the sky and the geometrical model of the source. We emphasize that transition to saturation is a 
gradual process. It is seen from Eq. (3.2.1) that to pass from an unsaturated regime { J/Jg ■C 1) to a 
heavily saturated one {J/Js » 1), the intensity needs to increase by about two orders of magnitude. 
Thus the saturation criterion J > Jg is no more than an order of magnitude criterion. 

There is one case for which we can simplify Eq. (3.2.9) without losing its exactness: the case where 
maser radiation is generated at the "optimum" density, i.e., the density that provides the maximum 
(unsaturated) value of \knet\ for the line in question. As was discussed in Section 3.1, the optimum 
density for maser amplification is much higher than the density of maximum inversion, and, as Figs. 5 
and 6 show, this optimum density is actually rather close to the density of thermalization. At such high 
densities, At in Eq. (3.2.9) can always be ignored as compared with (C( + C21). For 5 < n < 35 and 
Te = 10^ K, we use the following approximation: 

S = {Ct +C2i)^S- 10-^ n^Ne, (3.2.10) 

which is a simplified version based on Gee's etal. (1976) analytical approximations for collision rates 

and which is exact to ±5%. Substituting for the approximation (3.1.6) for A^™°^, and taking A21 « 
6.3 • 10^ nj"^ (s~^) and ~ 6.6 • 10^^ n^^ (s~^) — common approximations valid for a-lines with n 1 
— and finally neglecting At in Eq. (3.2.9) as compared with (Ct + C21), the equation reduces to 



The previous assumption, that the observed masing lines arise in regions of "optimum" density, is 

of course an assumption which can be wrong. One reason may be the interaction of the adjacent masing 
transitions due to saturation. Suppose some line 2^1 mases at its optimum density and reaches strong 
saturation. We can ask how such saturation might affect the trends of masing in the next transition 
higher, 3^2. If the 2^1 line were not saturated, the 3^2 line would be if not thcrmalizcd, at this 
high density, at least amplifying less strongly than its own optimum (lower) density would predict (see 
Fig. 5 for an illustration). According to Eq. (3.2.1), however, saturation of the 2^1 transition means 
that A7V21 decreases due to the transfer of some population from level 2 to level 1, and implies that the 
population of the lower level of the 3^2 transition decreases. In other words, saturation of the transition 
2^1 should increase the inversion of the transition 3 — > 2. In thermodynamic language, saturation of 
the transition 2^1 opens a new channel of sink to the maser on the transition 3^2. thereby increasing 
the efficiency of its pumping. As a result, the 3^2 line will be amplified to a higher intensity. If it in 
turn saturates, it will increase in the same manner the inversion of the next higher transition, etc. The 
efficiency of this upward diffusion of the rcenforced piimp sink decreases, however, at every step, until 
eventually, at some higher value of n, the steeply increasing relative rate of thermalizing collisions will 
remove the population inversion. 

Thus, a saturated high frequency masing line can "attract" an adjacent group of lower frequency 
lines to mase in a region of density higher than their unsaturated "optimum" densities. We will apply 
this mechanism, in Paper 3, to interpreting some observed properties of the multi-line hydrogen maser in 
MWC349. In any particular case, however, this qualitative idea requires a more exact mimerical check — 
a calculation of level populations including induced radiative transitions in the lines which are supposed 
to be saturated masers. 



4. HYDROGEN "DASARS" 

The effect of overcooling, demonstrated for hydrogen Rydberg levels in Fig. 4a and 5 and discussed 
in Section 2.3, is actually a well known phenomenon for radio lines of some interstellar molecules (H2CO, 
OH) . Observationally it appears as an enhancement of a line seen in absorption. The observational effect of 
overcooling was called "DAS AR" ( "Darkness Amplification by the Stimulated Absorption of Radiation" ) 
by C. Townes and A. Shawlow. Overcooling is sometimes called "anti-inversion," and the corresponding 
enhanced absorption — "anti-maser." The analysis of non-LTE states of a quantum transition given in 
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Section 2.2 shows, however, that overcooling is really more an analogue of overheating than of inversion 
(see more in Strelnitski, 1983). 

We first show that for an homogeneous, ionized cloud with no external continuum background, the 
observation of a line in absorption, i.e., r < 0, requires /3i2 > 1 - overcooling of the corresponding 
transition. 

Prom Eq. (2.1.6), r < corresponds to 

V < , ^-7(---) . . (4.1.1) 
1 -exp(-Tc -Ti) 

It is easily seen that the right-hand side of Eq. (4.1.1) is always < 1, whatever the value (or sign) of n, 
and hence rj < 1, which is only possible according to Eq. (2.1.2) if /3i2 > 1. 

Thus, in the framework of this simple model, the very observation of a recombination line in absorp- 
tion on the proper free- free continuum of the cloud is a direct indication that the transition is overcooled. 
More complicated circumstances, however a hotter external continuum background for example, or 
some peculiar morphology of the cloud, like perhaps if the kinetic temperature decreases radially out- 
ward — can produce an absorption recombination line without overcooling. Usually such cases can be 
identified by a detailed analysis of the source. 

Unfortunately, the reverse statement, namely that P12 > 1 (overcooling) always results in r < (a 
line seen in absorption), is not true. To see when overcooling does produce an absorption line in our 
model, we consider three particular cases: 

(1) Tc < 1, t; < 1; 

(2) Te<l, t;>1; 

(3) Tc>l. 

In the first case, expanding the exponents in Eq. (2.1.6) to first order terms and using Eqs. (2.1.2) 
and (2.1.4), we have: 

r^rjlf-l = ^b,. (4.1.2) 

Since all the terms on the right-hand side are positive, r can only be positive in this case. In other 
words, only emission recombination lines can be observed from an optically thin cloud with no external 
continuum background. 

Case (2) can occur for high frequency recombination lines (n<30), for which k'^ <C fc; (see Fig. 10). 
Ignoring k[. in Eq. (2.1.2) and expanding the denominator on the right-hand side of Eq. (2.1.6) to first 
order, we have, as the condition of r < for this case: 

To > 1^^^ . (4.1.3) 

pi 2 

Thus, for high frequency lines overcooling can produce an absorption line even in a cloud optically thin 
in the continuum, however, the line optical depth of the cloud needs to be ^1 and its continuum optical 
depth should satisfy the condition above (4.1.3). 

Finally, we consider case (3), first in its asymptotic form, when Tc » 1. Eq. (2.1.6) reduces then to 
r w 7/ — 1, and the condition r < becomes 77 < 1 which, according to Eq. (2.1.2), is identically satisfied 
if /3i2 > 1. Thus, for a homogeneous cloud optically thick in continuum an overcooled recombination line 
will always be observed in absorption (whatever the value of t; ) . 

Now we show that this is true even for a moderately optically thick cloud (tc^ 1). In this case both 
numerator and denominator in Eq. (2.1.6) are less than 1 and greater than (1 — e""'^) « 0.63. Thus 
the condition r < in this case becomes r] < 0.63. After some algebra on equation (2.1.2) for -q, and 
remembering that 62/&1 < 1 when /3 > 1, the r < condition becomes: 



k'lbi > 0.6 



k'c ~ (3i2 - 1.6 ■ 



(4.1.4) 



Taking /3i2 (at its maximum) from Fig. 5, bi from Fig. 2, and fc;,^^ from Fig. 10, it is easy to see that 
Eq. (4.1.4) is satisfied for all the transitions which can be overcooled in principle, n^;20 up to at least 
n — 56, seen in Fig. 5. Thus, for all these transitions (and in fact for higher n transitions as well) 
overcooled lines will be observed in absorption, if Tc J 1. 
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Observation of an absorption line; on the proper free-free continuum of an HII region is the most 
direct way to reveal overcooling, but it is not the only possible way. At low optical depths (Case 1) 
there is another obvious way: to observe more than one line in the interval of n where overcooling is 
theoretically predicted and to determine the decrement of r with respect to n. According to Eq. (4.1.2), 
in the optically thin case r oc 6„, hence, the observed sign of the decrement of r should show whether or 
not overcooling takes place. 

5. PROSPECTS 

In this Section we address the observational prospects for HRL masers, lasers, and dasars. Our anal- 
ysis of masing shows that there are at least two preliminary questions to be addressed when considering 
an object, or a class of objects, as a potential candidate: (1) what are the densities of ionized hydrogen in 
the source? and (2) what is the linear scale, or, more exactly, the column density per thermal line width, 
at each density? The answer to the first question shows us which lines might be inverted or overcooled in 
the source, the answer to the second — whether or not these non-LTE lines can become high gain masers 
or detectable dasars. 

We examine now, with these two test questions, several classes of objects. This list of candidates 
should not be regarded as a complete one: we only give several examples. We do not discuss here the 
objects similar to MWC349 — hot emission line stars with neutral disks ionized on their surface. Some 
observational prospects regarding this class of sources will be presented in Paper 3. 

5.1. Compact HII Regions 

Simple estimates show that among the compact HII regions associated with the sites of massive star 
formation, only the densest , "ultracompact," ones can have optical depths in HRL close or surpassing 
unity by absolute value as required for masing. Typical electron densities in ultracompact HII regions 
are ~ 10**-10^ cm~^, but in some of them an average density is estimated to be as much as Ng w 
3 • lO^cm"^ (Wood and Churchwell 1989; Kurtz etal. 1993). Figure 8 shows that the lines in a broad 
interval, from w 35q: to « 60q:, have comparably high |A:„ei| at this density, which may reach 10^"^ cm""'^, 
if the temperature is somewhat lower than 10,000 K. Since many of these ultracompact HII regions have 
dimensions L^IO^^ cm, the maser gain \knet\L may surpass unity in these lines, even if the velocity 
dispersion reduces the effective value of \knei\ by a factor of 1.5 2. Thus, some increase of intensity and 
line narrowing due to masing is anticipated in this interval of a-lines from the most compact HII regions. 

Some of these regions are highly inhomogeneous. Since chances of high-gain masing increase, in 
general, with the increasing density, masing in lower n lines is also quite possible from local density 
enhancements, even if such zones of higher density are considerably smaller than the dimensions corre- 
sponding to the average density of the object. 

HII regions associated with evolved objects, for example planetary nebulae, have in general lower 
densities. Nevertheless studies of forbidden line transitions (e.g., Osterbrock, 1989) show that in some 
cases like IC4997 there are regions with densities up to 10^ cm~^. 

Effects of masing can be studied by comparing the intensities and widths of several a-lines over a 
broad range of n, and also by comparison with the /3-lines from nearby n (see Paper 3). The analysis 
would help specify more exactly the physical parameters of the HII region. 

5.2. Wolf-Rayet Stars 

The idea to search for HRL line masers in Wolf-Rayet (WR) stars originated with H.E. Matthews 
(1992; private communication). These stars possess strong winds and expanding ionized circumstellar 
disks which might promote population inversions and maser amplification similar to those in the disk 
of MWC349. However, the first attempts to find mm and submm hydrogen masers in WR stars gave 
negative results (Matthews 1994; private communication). 
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The most probable reason is insufficient optical depth in the lines (SOa - 26a) studied. They require 
A^e ~ 10''-10* cm"'^ (see Table 1) for effective masing — but these densities are not sustained in large 
enough column densities in a WR disk outflow. If we approximate the ionized envelope of a WR star as 

a spherically symmetric, > 1000 kms^^ constant velocity (thus N,, oc outflow with a mass loss rate 

of < lO"** M0 yr^^, it is easy to show that the linear scale of the region of density Ng will be R{Ne) < 2 • 

lO^'^A'^e cm. At the optimum density of the H26a line, for example, where iV™"^ « 1 • 10* cm~^, 

we find i?(7VJ"°^) < 2 • 10^^ cm. Table 1 gives fc™°^ 10"" cm"! for this fine, but, in calculating 

the optical thickness we must reduce this value by the ratio of the radial velocity dispersion in a WR 
envelope, AV ^ 10'^ kms~^, to the thermal line width, vth ^ 20kms^^ (see Sect. 3.1). As a result, the 
net optical depth in the line reduces to T„et ~ knei -R(A^""'' ) (vih/AV) ^ — 10~^. In the same way we can 
show that |r„ei| ^ 1 for other mm and submm lines. Though |r„e/,| does increase with decreasing r?., it 
becomes ^ 1 only at 10. Thus, high gain lasers in hydrogen recombination a-lines from the expanding 
shells of WR stars are most promising for the IR lines with n < 10. 

At this high a frequency a different source of line emission can be even more effective at producing 
laser lines in WR stars: the atmosphere. In the standard model of a WR star the non-expanding 
atmosphere consists of a turbulent layer where Ne ~ 10^^-10^"^ cm""^ and which is ~ 5 • 10^^ cm thick. 
Such densities are optimum for lasing in lines (see Table 1). Multiplying the width of this layer 
by our calculated values of fc™"^, we find very large optical depths; with Ne ^ 10^^ cm^'', for example, 
the optical depth of Ba should be ~ —200! Even assiiming the lower limit to the density in this layer 
(deduced by Aller from optical spectroscopy) of Ai'e ~ 10^^ cm~^, we still have optical depths of —0.2, 
—3, and —2 for the Ba, H6a, and H7a lines respectively, with Ba being the most promising line if the 
density is higher than this lower limit. 

The observability of an HRL laser is, however, strongly contingent upon its geometry. If the solid 
angle of lasing is not small enough, this emission tends to be lost on the stronger, spontaneous emission 
background. This problem is discussed in more detail in Paper 3. 

5.3. Be Stars 

Smith et al. (1979) suggested that Be stars might be likely candidates for lasing in the infrared 
hydrogen lines, although at that time the lack of theoretical high density population data precluded a 
quantitative discussion. These complex objects have ionized Keplerian disks (Slettebak etal. 1992), high 
mass loss rates, photometric variability, and, most important for our considerations, inferred densities 
in the region emitting Balmer lines of Ne ~ 10^^ lO^'^cm^'', over dimensions ~ 10^^-10^^ cm. These 
parameters are similar to those in the turbulent layer of a WR atmosphere, and similar considerations 
hold here too, suggesting they are promising candidates for lasing in the near IR recombination a-lines 
(n<10). 

5.4. Active Galactic Nuclei and Starburst Galaxies 

Enormous masses of ionized hydrogen are present in the central parts of galaxies with active nuclei 
(AGN), while starburst galaxies, which may sometimes be the result of two galaxies in collision, often 
have bright compact HII regions outside the nucleus. In the case of AGN's we consider the possibilities of 
maser action in two types of regions (e.g., Osterbrock 1989). In the so-called "Broad Line Regions" (BLR) 
Ne ^ 10^-10^° cm""^, and may range between L ^ 0.01 0.1 pc (Seifert 1 galaxies) to ~ 1 pc (quasars). In 
the "Narrow Line Regions" (NLR) A^e ~ 10^-10^ cm~^, L ^ 1000 pc. However, the densities above do not 
apply to these linear scales but ruther to much smaller "clumps" or cloudlets — the filling factor should 
be much less than one. 

At the densities typical for BLR, Ha-lines with n « 10 — 25 have maximum local maser gain (Fig. 8). 
Taking the values of fc™""^ for these lines from Fig. 6 or Table 1, we see that linear scales of only ~ 10 AU 
(25a) to ~ 0.01 AU (10a) are needed to make these lines optically thick, if line broadening is thermal. 
If the individual cloudlets in BLR are of these dimensions or bigger, then laser emission in these submm 
and IR lines is anticipated, and would be seen as relatively strong narrow features centered at the radial 
velocities of individual cloudlets. It seems quite possible. For example, in the Scoville and Norman 
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(1988) model of the broad emission line clouds as photoionized mass-loss envelopes of giant stars, the 
Stromgrem sphere reaches down to a radius of ~ 10^'* cm in the stellar envelope. With the density there 
being ~ 10^ cm~^, the 10a-15a lines should have high unsaturated laser gains — on the order of tens. 

The NLR gas seems to be a less likely prospect for high-gain masing in HRLs. The densities of the 
NLR gas are optimum for masing in high n a-lines — n>50. Densities ~ lO^'cm"'^ have linear scales 
up to ~ Ipc in AGNs (viz., the ionized accretion disk "halo" postulated by Scoville and Norman, 1995). 
According to Fig. 8, the 50a - 80a lines have the highest \knet\ at this density, but it is easy to see that 
maser gains of only <0.1 are achievable. The regions with A^e ~ lO^cm"'^ can be much larger, up to 
Ikpc, but considering the velocity dispersion in such a large region, we have to reduce the optical 
depth at each frequency by a factor of AV/vth ~ 50, and Tnet turns then out to be < 1 in this case too. 

Maser effects in mm HRLs are anticipated in the star burst galaxies whose archetype, M82, was in 
fact the first external galaxy where a radio recombination line (H166a) had been detected (Shaver etal. 
1977). Kronberg etal. (1985) postulate dense (Ag ^ lO^cm"^) ionized shells surrounding supernovae 
which had exploded inside dense molecular clouds of the M82-like starburst galaxies. According to 
Fig. 8, a-lines with n « 20-35 have the highest maser gains at this density. With the typical linear scale 
for the shells being ^ 3 • 10^^ cm. Fig. 8 shows that the gain can surpass unity, if large-scale motions do 
not broaden the HRL too much, as compared with the thermal width, and if Tg is not much higher than 
10,000 K. Though the both latter assumptions can be wrong, some indications of narrow, presumably 
weakly masing components in the 36a line from M82 has recently been detected by Smith etal. (1995). 

5.4- Hydrogen Dasars 

Prom an observational point of view the major difficulty with the overcooling phenomenon in HRLs, 

as compared with inversion, is the low upper limit for Ag: overcooling disappears at Ne^lO^ cm~^ (see 
Sect. 4). At this limiting density hydrogen a-lines with n « 20-25 are overcooled, and the Case 2 (cf., 
Sect. 4) applies. Since /3i2 is only slightly greater than 1 for these lines, Eq. (4.1.3) shows that n must be 
> 1 to observe an absorption line on the proper continuum of an HII region. The typical value of fc; for 
these lines, at Ng ~ lO^cm"^, is ki ~ 10~^°cm~^ (Fig. 10), and thus very large Hll region dimensions, 
of > 300 pc, are needed for t; > 1. HH regions with such parameters are unknown. 

At lower densities, which might make it easier to find larger regions, the maximum of overcooling 
shifts to the higher n lines. The values of 6„, /3i2 and the oscillator strengths of the lines all increase, but 
the Ag factor decreases faster, so the required linear scale increases. Of known objects, we believe only 
AGN's can be considered today as possible candidates for direct observations of hydrogen dasars. If the 
regions of intermediate density, Ag ~ lO^cm"^, have dimensions comparable with NLR's (up to ~ 1000 
pc), a-lines from n « 30-35 may acquire 77 ~ 1 and appear as dasars. 

The most realistic way to detect the effect of overcooling in hydrogen, however, is to observe more 
than one line in the n range where overcooling is predicted at a given density, and to determine the 
decrement of r with n, as described in Sect. 4. This way of inferring overcooling is applicable to any 
galactic HH region if the observations are done in the proper frequency range. 

6. CONGLUDING REMARKS 

Hydrogen, the simplest and the most abundant atom in the Universe, is also the first known atomic 
astrophysical maser. Perhaps the most striking feature of this maser, as compared with molecular astro- 
physical masers, is its obstinate uniqueness: six years after its discovery in MWC349 there is still only one 
known strong hydrogen maser source, despite intense attempts to find another. In contrast, practically 
every newly discovered molecular maser line was quickly detected in a multitude of similar sources. 

Any useful theory of the hydrogen maser phenomenon must directly address this striking fact. 

In this paper we have attempted to show that, although hydrogen may be common, the physical 
conditions necessary for creating a high gain HRL line maser are rather special: for every line (or narrow 
group of lines) only a limited range of density is able to producing the high negative values of the local 
absorption coefficients required. Furthermore the region with this special density must be adequately 
extended, and homogeneous enough in its radial velocity, to create an optical depth < — 1. The velocity 
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dispersion parameter is an especially important one. One of the advantages of molecular masers is that 
they are located in relatively cold, neutral gas clouds and are produced by relatively heavy particles, so 
the width of the local thermal velocity profile is, in general, only ~ lkms~^. Together with the high 
density of these cold clouds, such a small velocity dispersion can produce over even a short distance a 
high negative optical depth at the line center, and thus very strong masing. A HRL maser, on the other 
hand, is generated in a hot ionized gas and is radiated by the lightest particle, whose thermal velocity 
dispersion is therefore very high, ~ 20kms~^. As a result the local absorption coefficient is relatively 
low, regardless of the degree of inversion, and so long distances are required to create r < — 1. But 
longer distances bring an additional velocity dispersion to the picture due to their large scale motions — 
expansion, contraction, and the like, which broaden the frequency profile of the amplified radiation even 
more, thus hindering strong masing. 

In such situations chance may become a decisive factor for creating an observable maser. In the case 
of MWC349, there arc at least two lucky coincidences: (1) a very high UV- luminous a star, perhaps 
in a particular outflow phase of its life, surrounded by a dense, massive neutral disk, resulting in an 
unusually dense, extended HII region on the disk's surface; and (2), a fortuitous projection of the disk to 
our line-of-sight as edge-on. A Keplerian velocity field seen edge-on is a better "organizer" of the atoms' 
radial velocities than an outflow — that is why we observe stronger mm and submm masers from the 
disk of MWC349 than from its outflow (cf. Paper 3). 

Whatever the particulars of HRL line masers and lasers, nature is fabulously inventive. Conditions 
favorable to the creation of rather strong hydrogen masers and lasers, or even abnormally strong ab- 
sorbers (dasars), might occur in several classes of astrophysical objects, both galactic and extra-galactic. 
The strong hydrogen maser in MWC349 was discovered accidentally, during a routine search for more 
recombination lines in this well-known emission line star, at a time when the paradigm of "decrease of 
stimulated emission effects toward higher frequencies" was strongly held. Now, thanks to the example of 
MWC349 coupled with a thorough analysis of high density hydrogen level populations, wc can be much 
more optimistic about finding this intrinsically interesting and informative phenomenon elsewhere in the 
Universe. 
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FIGURE CAPTIONS 



Fig. 1. To the classification of non-LTE states of a quantum transition. 

Fig. 2. Solid lines: 6„ (n) at Tg = 10** K and at different electron number densities (cm~^), 
indicated near the curves. Numerical data from Storey and Hummer (1995). Broken line: 6„ (n) at 
Ne = ("radiative limit"), calculated with the program of one of us (V.P.) mentioned in text. 

Fig. 3. To the explanation of the overcooling mechanism: bni (l) for n = 30 — 45, at Tg = 10^ K 
and A^e = lO^cm"^. Note the rearrangement of the 6„; factors in favor of the high-? sublcvcls with 
the decrease of n, due to the increasing unblurring of the ^-structure. Numerical data from Storey and 
Hummer (1995). 

Fig. 4. a — /3n,n+i(n) at Te = 10'' K and at A^e = 10"^ — lO'' cm~^ (numbers near the curves). 
b — The same for A^e = 10^ — 10^^ cm~^ . Numerical data from Storey and Hummer (1995). 

Fig. 5. Pn,n+i{Ne) at Tg = IO^'K for different na lines (n is indicated near the curves). Numerical 
data from Storey and Hummer (1995). The bars, in the lower part of the figure, show the location and 
the width (FWHM) of the fc; {Ng) and knet (-^e) profiles for the 36a line, according to our calculations 
(compare with Fig. 6b). 

Fig. 6. a-i — knet (-^e) for some a-lines at = 10* K. j — knet i^e) at Te = 10* K for some /J-lines. 

Fig. 7. Solide lines: log |fc„et| (log A^e) for higher n a-lines (n is indicated near the curves), for 
Te = 10* K. Broken lines loci of the maximal |fc„et| for individual a-lincs, for Te = 10* and 5 • 10^ K. 
Dotted line - the envelope of the \knet \ (-^e) plots for Te = 5 • 10*^ K (the plots are not shown). 

Fig. 8 Log \knet \ {'n) for different electron densities, for Te = 10* K (solide lines) and 5 • 10^ K (broken 
lines). 

Fig. 9. Dots: Graphical representation of the A^™°^ (n) dependence, with the numerical data from 
Table 1 for Tg = 10* K. Dashed line: approximation (3.1.2) for 5 < n < 40. 

Fig. 10. The LTE line and continuum absorption coefficients, k'l and k'^, as functions of n, for 
Te = 10* and 5 • 10^ K. 



TABLE CAPTION 

Table 1. Optimum densities of unsaturated mascr amplification, N^"^^ and the values of A;™"^ for 
hydrogen recombination a- and /3-lines at Te = 10,000 and 5,000K. 
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